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Utilizing a silicon nanobeam optomechanical crystal, we investigate the attractor diagram arising from the 
radiation pressure interaction between a localized optical cavity at X c = 1552 nm and a mechanical resonance 
at co m /27l = 3.72 GHz. At a temperature of 7 1 « 10 K, highly nonlinear driving of mechanical motion is 
observed via continuous wave optical pumping. Introduction of a time-dependent (modulated) optical pump is 
used to steer the system towards an otherwise inaccessible dynamically stable attractor in which mechanical self¬ 
oscillation occurs for an optical pump red-detuned from the cavity resonance. An analytical model incorporating 
thermo-optic effects due to optical absorption heating is developed, and found to accurately predict the measured 
device behavior. 


The field of optomechanics, concerned with the interaction 
of an optical cavity and a mechanical resonator |T|, has been 
of recent interest for its promise for use in sensors EH, non¬ 
linear optics nm, and demonstrations of macroscopic quan¬ 
tum mechanics nm. To lowest order, the mechanical dis¬ 
placement linearly modulates the frequency of the optical res¬ 
onance in a cavity-optomechanical system. This, however, 
gives rise to an inherently nonlinear phase modulation, and 
through radiation-pressure backaction on the mechanical el¬ 
ement, yields nonlinear system dynamics (8)- Much of the 
previous work has focused on the linearized regime where the 
interaction with the optical field still gives rise to a host of in¬ 
teresting phenomena such as a modified spring constant HE 
damping or amplification of the mechanics iflOi . and EIT-like 
slow-light effects [11,12 ]. Recently, several experiments have 
pushed into the quantum regime using back-action cooling in 
the linearized regime to damp nanomechanical resonators to 
near their quantum ground state of motion [ 13, .14]. 

In this work, we instead demonstrate new features and 
tools in the nonlinear regime of large mechanical oscilla¬ 
tion amplitude. Previous experimental works have shown 
that a blue-detuned laser drive can lead to stable mechanical 
self-oscillations [ 15-18], or even chaotic motion fl9l . The¬ 
oretical predictions of an intricate multistable attractor dia¬ 
gram (8) have so far eluded experimental observation, except 
for the elementary demonstration of dynamical bistability in a 
photothermally driven system [201. In the present work, we 
are able to verify the predicted attractor diagram and further 
utilize a modulated laser drive to steer the system into an iso¬ 
lated high-amplitude attractor. This introduces pulsed control 
of nonlinear dynamics in optomechanical systems dominated 
by radiation pressure backaction, in analogy to what has been 
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shown recently for a system with an intrinsic mechanical bista- 
bility EQ. 

We employ a one-dimensional (ID) optomechanical crys¬ 
tal (OMC) designed to have strongly interacting optical and 
mechanical resonances E2 The OMC structure is created 
from a free-standing silicon beam by etching into it a peri¬ 
odic array of holes which act as Bragg mirrors for both acous¬ 
tic and optical waves [23]]. A scanning electron micrograph 
(SEM) of an OMC cavity is shown in Fig.[l^ along with finite- 
element-method (FEM) simulations of the co-localized opti¬ 
cal (Fig.[]J>) and mechanical (Fig. |TJ:) resonances. To reduce 
radiation of the mechanical energy into the bulk, the OMC is 
surrounded by a periodic ‘cross’ structure which has a full 
acoustic bandgap around the mechanical frequency (Fig. [T^, 
green overlay) ff24l . 

The experimental setup is shown schematically in Fig. [TJi. 
The silicon chip containing the device is placed into a helium 
flow cryostat where it rests on a cold finger at T ^4K (the 
device temperature is measured to be 7J, « 10 K). Input laser 
light is sent into the device via a tapered optical fiber, which, 
when placed in the near-field of the device, evanescently cou¬ 
ples to the optical resonance of the OMC [25]. The trans¬ 
mitted light is detected on a high-bandwidth photodiode (Dl) 
connected to a real-time spectrum analyzer (RSA). We also 
employ an electro-optic modulator (EOM) in the laser’s path 
to resonantly drive the mechanical resonator. Finally, we can 
send in a low-power, counter-propagating probe laser whose 
detected spectrum (D2) is used to measure the mechanical am¬ 
plitude and the pump-cavity detuning. Using this set-up the 
optical resonance of the device studied in this work is mea¬ 
sured to be at at X c = 1542 nm, with intrinsic (taper loaded) 
energy decay rate of k*/ In = 580 MHz (k/271 =1.7 GHz). 
The breathing mechanical mode, shown in Fig. [Tjc, is found 
to be at co m /27t = 3.72 GHz with bare energy damped rate of 
Y//27U = 24 kHz. 

The coupling of the optical resonance frequency to the 
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FIG. 1: (a) SEM of the optomechanical nanobeam surrounded by 
phononic shield (green), (b) FEM-simulated electromagnetic energy 
density of first-order optical mode; white outline denotes edges of 
photonic crystal, (c) FEM-simulated mechanical mode profile (dis¬ 
placement exaggerated). In (b) the colorscale bar indicates large (red) 
and small (blue) energy density, whereas in (c) the scale bar indicates 
large (red) and small (blue) displacement amplitude, (d) Simplified 
schematic of experimental setup. WM: wavemeter, A<|): electro-optic 
phase modulator, OMC: optomechanical crystal, Dl: pump light de¬ 
tector, D2: probe detector, VOA: variable optical attenuator, PM: 
power meter 


mechanical displacement yields the interaction Hamiltonian, 
//int = hgoa^dx where a (. x ) is the optical (mechanical) field 
amplitude, and go is the vacuum optomechanical coupling rate. 
The physical mechanical displacement expectation is given by 
x = Xzpf(x), where the zero-point amplitude of the resonator is 
x Z pf = (^/2m e ffC0 m ) 1 ^ 2 = 2.7 fm (estimated using a motional 
mass m e ff = 311 fg calculated from FEM simulation). Uti¬ 
lizing a calibration of the per-photon cooling power ED we 
find that go/2ft = 941 kHz. These device parameters put our 
system well into the sideband resolved regime K/co m <C 1. 

The classical nonlinear equations of motion for the mechan¬ 
ical displacement (jc) and the optical cavity amplitude (a = (a)) 
are, 


x(t) = -y ix(t) - c 0 2 m x(t) + 2® m g 0 x zp f\a(t)\ 2 , (1) 


a(t) = 


— ^ + z' ( Al + -^—x(t ) 

■^zpf 


a (t) + \j y a m , (2) 


where a m = y.P m //i<»/ is the effective drive amplitude of the 
pump laser (input power P[ n and frequency Cl>l), K e /2 is the 
fiber taper input coupling rate, (D c is the optical cavity reso¬ 


nance frequency, and = co^ — C0 c . Since we are interested in 
the regime of self-sustained oscillations, where the motion of 
the oscillator is coherent on time scales much longer than the 
cavity lifetime, we can take the mechanical motion to be sinu¬ 
soidal with amplitude A, x(t) = Asinco m E The optical cavity 
field is then given by, 



where <F(7) = — p m cosco m £ and a n = 

J n (Pm) / (k/2 + i OcD m — Al)). Here J n is the Bessel function 
of the first kind, ft-th order, and its argument is the unitless 
modulation strength p m = (Ago) / (*zpftt>m)- For Pm ^ 1 only 
the terms oscillating at the mechanical frequency, co m , are 
appreciable, so the interaction can be linearized, and only the 
first-order radiation pressure terms are present. However, for 
P > 1 the higher harmonic terms at each n co m have significant 
amplitude and backaction force. 

The thermal amplitude is too small to enter the nonlinear 
regime in our devices (Pth ~ 0.01), however, backaction from 
the pump laser can provide amplification to drive the mechan¬ 
ical resonator into the high-P, nonlinear regime. The result¬ 
ing mechanical gain spectrum in the amplitude-detuning plane 
(the attractor diagram) can be solved for by calculating the en¬ 
ergy lost in one mechanical cycle (Pf r i c = m e ffy \ (x 2 )) and com¬ 
paring it to that gained (or lost) from the optical radiation force 
(F ra d = (%oAzpf) ( \a\ 2 x )) 0 . Figure Ua shows a plot of the 
gain spectrum for the parameters of the device studied here 
with a laser pump power of P m =151 juW. Imposing energy 
conservation, Prad/Ffric = +1, yields the steady-state solution 
contour lines. Although the entire contour is a physical solu¬ 
tion, the equilibrium is only stable when the power ratio de¬ 
creases upon increasing the mechanical amplitude, Jg ^ < 0 
(i.e. stability is found at the ’tops’ of the contours) (8). 

In the device studied here there is a thermo-optic frequency 
shift of the optical cavity caused by heating due to intra¬ 
cavity optical absorption. The thermal time constant of the 
device structure is slow relative to the optical cavity coupling 
rate, but fast compared to the laser scan speed. Absorption 
heating can thus be modeled as a shift of the laser detuning 
proportional to the average intra-cavity photon number (fi a ) 9 
A l = Al,o + Ct 0 n a , where Al,o is the bare laser-cavity detuning 
in absence of thermo-optic effects. The per photon thermo¬ 
optic frequency shift of the optical cavity is measured to be 
Cto/271 = —216 kHz. Including this effect, the shifted con¬ 
tours are shown in Fig. |2j) as a function of the bare detun¬ 
ing A^o- The solid lines with arrows indicate the expected 
path traversed by the mechanical resonator during a slow laser 
scan from lower to higher laser frequency (left to right) at each 
power. The dashed lines are contours that are either unstable, 
or unreachable by this adiabatic laser sweep. 

We first explore the lowest-lying contour of the attractor 
diagram by measuring the optical transmission as the pump 
laser is tuned from red to blue across the optical cavity res- 
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FIG. 2: (a) Calculated gain spectrum for the OMC in the amplitude-detuning plane. Color scale indicates the ratio of power input to that lost 
from friction (P r ad/^fric — 1) at Pi n = 151 jliW. Positive values are regions of mechanical self-oscillation. Solid line curves indicate power- 
conserving solution contours at selected input powers: 0.65 juW (white), 6.5 juW (grey), 151 juW (black), (b) Same as (a) with contours 
now shifted by estimated thermo-optic effects. Solid line curves indicate the path taken by the mechanical oscillator during the laser sweep. 
Dashed lines are contours which are either unstable or unreachable by a slow sweep of laser detuning, (c) Image plot of the measured optical 
transmission spectrum versus laser detuning and power, (d) Image plot of the theoretically calculated transmission spectra including thermo¬ 
optic shifts and a slow drift in the optical resonance frequency over the course of the measurement from low to high power. Spectra in (c) and 
(d) are scaled at each power level to span the range 0—1. (e) Plot of the normalized optical transmission from scans in (c) at Pi n = 0.12 juW 
(top), 0.65 fji W (center), 151 juW (bottom). Blue points are measured data and red curve is the theoretical model, (f) Power spectral density of 
detected signal near the mechanical frequency for P[ n = 151 ju W, showing frequency shifts of the mechanical mode, A m = co m — co m5 o, from its 
bare frequency 0d m ^/2n = 3.72 GHz. Color scale is detected power density in dBm/Hz. (g) Total integrated power of spectra in (f). Measured 
data are plotted as green circles, with the theoretical model (up to a scale factor) shown as a solid red curve. The red arrow in each plot 
indicates the laser scan direction. 


onance with a fixed optical input power. A dip in transmis¬ 
sion indicates that light is entering the cavity and being lost 
through absorption or scattering. At low optical input powers 
(Pin < 0.3 /iW), only a single resonance dip associated with 
the bare optical cavity is observed. Upon increasing the laser 
power, radiation pressure backaction amplifies the thermal mo¬ 
tion of the mechanical resonator beyond threshold and into a 
large amplitude state. When this occurs a large fraction of the 
intra-cavity photons are scattered, resulting in additional trans¬ 
mission dips near each detuning A^o = where threshold 
is reached. Physically, mechanical oscillations at the n -th side¬ 
band detuning are generated by a multi-photon gain process in¬ 
volving n photon-phonon scattering events. This stair-step be¬ 
havior is seen in the measured transmission spectrum of both 
Rg# and Fig. ET The theoretically calculated spectra in¬ 
cluding thermo-optic effects are shown in Fig. [2ji, in excellent 
agreement with the measured spectra after taking into account 
a slow drift in the cavity resonance frequency as the measure¬ 
ments were taken from low to high power. 

Figure (2]f shows the radio-frequency noise power spectrum 
near the mechanical frequency of the optical transmission pho¬ 
tocurrent at the highest measured input power (Pin = 151 \i W). 
We note that backaction effects blue-shift the mechanical res¬ 
onance frequency by an appreciable amount (^ 2 MHz) from 
its intrinsic value of C0 m /27C = 3.72 GHz. This frequency shift, 


A m , depends on the laser detuning in a more intricate fashion 
than a linear calculation of the optical spring effect would sug¬ 
gest. A measure of oscillation amplitude can be extracted from 
the total power in this mechanical sideband, and its measured 
dependence on detuning is shown in Fig. [2j*. The measured 
total transduced power oscillates due to the nonlinearity of the 
detection process, which is nicely captured by the theoretical 
model, obtained without fit parameters except for an overall 
scale factor (Fig. [2^, red curve). 

It is readily apparent from Fig. [2^ that at large optical pow¬ 
ers (black contour) there are a number of isolated attractor con¬ 
tours at higher oscillation amplitudes. Here we utilize exter¬ 
nal time-dependent driving of the mechanical mode to explore 
the lowest-lying isolated attractor on the red side of the opti¬ 
cal cavity (A l < 0), where the linearized theory predicts only 
damping of the mechanical mode. An electro-optic modulator 
(EOM) is utilized to phase-modulate the incoming light field, 
resulting in an oscillating force inside the cavity which drives 
the mechanical resonator towards higher amplitudes. The ex¬ 
perimental sequence is displayed in Fig. [3^. We start with the 
pump laser switched on at a power of P m = 43 ju W, the laser 
detuned to the red side of the cavity resonance, and the phase 
modulation off (P E0M = 0), which initializes the mechanical 
resonator into a cooled thermal state with p m « 0. The EOM 
phase modulation is then turned on which rings up the me- 
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FIG. 3: (a) Schematic showing the state preparation and measurement sequence used to explore the higher lying attractors, (b) Measured 
time-domain signal of one quadrature of the mechanical amplitude during the turn-on of the EOM drive (t = 0). For clarity, the signal is mixed 
down from the mechanical resonance frequency at 3.7 GHz to 150 kHz. (c) Measured mechanical amplitude (mixed down to 500 kHz) during 
the turn-off of EOM drive modulation at t = 0. Here we show a case where the system remained trapped in the high-amplitude state, (d) 
Transmission scans of the counter-propagating weak probe laser during each sequence of the measurement. Blue points are measured data and 
red curves are fits to the data using Eq.[3] A p = ti) p — co c . (e) Mechanical power spectrum (green points) near the mechanical frequency during 
each sequence of the measurement. In the left panel the resonator is in a cooled thermal state and the red curve is a fit used to extract mechanical 
resonator parameters, (f) Plot of the normalized gain spectrum in the detuning-amplitude plane with overlaid stable solution contour (black 
solid curve) at P m = 43 juW. Color scale is (P r ad/Pfric — !)• Dashed black curve indicates unstable portion of contour. Red data points indicate 
initial (p m , Al) and grey data points indicate the final values, for selected instances of the experimental sequence as shown in (a). White arrows 
connect initial/final pairs but do not indicate the actual path taken by the system. 


chanical resonator. Following ring-up, the pump laser is tuned 
to a starting detuning of A l, completing the initialization se¬ 
quence. Finally, the modulation is switched off (P E0M 0) 
and the system is allowed to relax into a final mechanical os¬ 
cillation amplitude and laser-cavity detuning. 

In order to determine p m and A L at each stage of the 
above procedure an additional counter-propagating weak op¬ 
tical probe beam of frequency co p is scanned across the cav¬ 
ity (see panels in Fig. [3ji). When the mechanical amplitude 
is large (p m > 1) the standard single resonance dip is trans¬ 
formed into a multi-featured spectrum, with resonance dips at 
the motional sidebands (A p = co p — (D c = ±^co m ). A fit to the 
probe spectrum is performed using Eq. [3] with p m and CD c as 
free parameters (Fig. Eh red curves). Repeating the measure¬ 
ment for different initial states and recording the resulting final 
states reveals the flow in the underlying attractor diagram. A 
representative subset of the results are plotted in Fig. [3]f. We 
find that for a narrow range of initial conditions, after the mod¬ 
ulation is switched off the system remains trapped close to the 
predicted top of the higher amplitude attractor at p m « 3.5 and 
Al/ co m « —0.7. For more negative initial detunings or lower 
initial mechanical amplitudes, the system relaxes into the triv¬ 
ial low-amplitude state or gets caught on the lowest-lying con¬ 
tour explored in Fig. [2] (for detunings A^/cOm > —0.5 the sys¬ 


tem could not be stably initialized due to the thermo-optic ef¬ 
fect). 

The results presented here represent an initial exploration 
of the nonlinear attractor diagram of an optomechanical sys¬ 
tem where the dominant nonlinearity is that of the radiation 
pressure interaction. Due to the limited drive amplitude of the 
electro-optic modulator used in this work (P E0M < 3.5), we are 
limited to exploring only the lowest red-side attractor. With 
the ability to apply larger drives, or to rapidly detune the laser, 
it should be possible to reach higher-lying islands, and more 
fully explore the attractor diagram shown in Fig. [2^. Further 
understanding of the latching effects in these measurements 
should also pave the way to exploiting them for use in metrol¬ 
ogy experiments as the dynamics that govern whether the os¬ 
cillator stably latches into an attractor can be a very sensitive 
function of the oscillator’s displacement (8), thus yielding a 
precise measurement of the oscillator’s environment or state. 
This latching also allows for systems with memory due to the 
hysteretic nature of the nonlinearity, as in tTR [2614291 . Finally, 
in future devices where the optomechanical coupling rate is 
larger, these same nonlinearities can lead to quantum mechan¬ 
ical effects which have thus far only been explored theoreti¬ 
cally 1301 . 
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